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Abstract 

We present the most general supersymmetric amplitude for K° — K° and 
mixing resulting from gluino box diagrams. We use this amplitude 
to place general constraints on the magnitude of flavor- changing squark 
mass mixings, and compare these constraints to theoretical predictions 
both in and beyond the Minimal Supersymmetric Standard Model. 
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Flavor changing neutral currents (FCNC's) provide an important test of the 
radiative structure of the Standard Model and a sensitive probe of new physics 
beyond the Standard Model. 1 Indeed, new particles proposed by alternative the- 
oretical frameworks can generate observable FCNC effects, even for new particle 
masses well beyond the range of present and proposed accelerators. For many such 
frameworks, K° — K° and B° — B° mixing provide the most sensitive experimental 
probe of these FCNC's. 1 

During the past eight years, there have been a number of theoretical papers 
exploring FCNC's in supersymmetry. 1-18 These papers generally disagree on the 
magnitude of supersymmetric FCNC's by as much as several orders of magnitude. 
These disagreements arise because there has been no systematic calculation of 
the relevant supersymmetric amplitudes for most FCNC processes comparable to 
those that have been performed in the Standard Model. 19 This is because the 
supersymmetric contributions are more complicated, involving many new particles 
(e.g., six complex charge 2/3 squark fields) which mix in a complicated way. 

In this paper, we present the exact Feynman amplitude for the dominant, 
gluino box contribution to K° — K and B° — B° mixing. This calculation may be 
organized in either of two frameworks: 20 first, using exact Feynman rules for mass- 
eigenstate squark fields with arbitrary off-generational squark-gluino couplings; 
second, as a perturbative expansion in small, off-generational squark mass inser- 
tions assuming diagonal squark-gluino couplings (and non-mass-eigenstate squark 
fields). The latter formulation affords the most convenient language both for ex- 
pressing experimental constraints and for comparing these constraints against dif- 
ferent theoretical predictions resulting from various supersymmetric models. 

We have verified that the result in the mass insertion framework can be system- 
atically derived from the result in the mass eigenstate basis by expanding around 
the universal squark mass. 20 This provides an important consistency check on the 
results obtained from both computational frameworks. Since, for Mq ^ 1 TeV, the 
magnitude of the expansion parameters, 5fh^ B /M~, consistent with experiment 
are all much less than one, an expansion to first order in these small parameters 
provides an excellent approximation. Apart from algebraic corrections, the main 
improvements over earlier computations 11 ' 12 are the clarification of three distinct 
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terms involving left-right mass insertions ( 15m 2 , 1 , \5rri 2 , 1 , 5m 2 , <5m| ) 

» » V L «£ s fl J L flfiSiJ ' d L s R d R s L ) 

and the identification of how the experimental constraints on, and theoretical pre- 
dictions for, the mass insertions scale with the supersymmetry breaking scale. 

The dominant supersymmetric contribution to AS = 2 processes results from 
four topologically distinct gluino box diagrams (Figures la - d). Neglecting exter- 
nal quark momenta compared to heavy internal squark and gluino masses, these 
diagrams gives rise to an effective interaction Lagrangian: 

a 2 { (brh\ \ 2 

+ 



with 
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f(x) = . 1 .- (-6 In a;- 18a;lnx-x 3 + 9x 2 + 9:r- 17) (2a) 
6(1 — x)° 

f(x) = - 1 J -^lnx-Qxlnx + lx* + 9x 2 - 9x - 1) (26) 
3(1 — x) b 

x = M? /M? (2c) 
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where Mq is the universal (or average) down-squark mass and the quantity [66f(x) + 
24xf(x)} — > — 1 for x — 1. The soft supersymmetry-breaking FCNC mass inser- 
tions Srhq (A, B = L, R) appearing in (1) are defined by 

C C -d\5m\ ASB S B (3) 

where the squark fields cLa-, §b reflect a "super KM" basis 20 in which the g, 7, Z 
gaugino couplings are flavor diagonal and the charged-current W mixing angles 
linking quarks and squarks are equal to the standard KM angles. 

One can easily verify that the interaction Lagrangian (1) reproduces the quark 
scattering amplitude in Figure 1, noting that the following operators give rise to 
associated amplitudes: 

[dalnPLS a ][dp^P L sp\ 

(4a) 

-> A = 2(d l3 (k 2 )-f^P L Sa(k 1 )) {d^ih^PLSsih)) (5 a p5 lS + 5 ai 5p 5 ) 

[diP A Sj][d k P B si] -> A = {dp(k2)PASa{ki)){d 1 {kz)PBS5{kA))8j a 5 i p8 kl 5i 5 

- (d^(h)P A s a (ki))(dp(k2)PBSs(h))Sja5i^5 k/ 35is 

(46) 

- (dp( y k 2 )PASs(h))(d> y (h)PBS a (ki))5js5ip5k~ / 5i a 
+ (d^(k 3 )P A s s (h)) (dp(k 2 )PBS a (k 1 ))5 j s5i 7 5 kl3 5i a 

The form of £^£=2 ^ n ^ ne " m i n i m al susy limit" where 5rh\ L dominates differs from 
previous analyses, 10 ' 11 which reported color octet operator structures in addition 
to the color singlet operator shown in (1). (A color octet left-left operator 

can always be converted to a color singlet operator by a Fierz transformation.) 
One can use the effective Lagrangian (1) to assess the effects of off-generational 
mass mixings (Jm^^of arbitrary chirality which are radiatively induced in minimal 
supersymmetry or which may appear in non-minimal supersymmetric models. The 
AB = 2 effective interaction Lagrangian can be read directly from (1) by the 
substitution 5m% — * 5m% , . 

The application of the low-energy effective AS = 2 (AB = 2) interaction 
Lagrangians derived above to K° — K° (B° — B°) mixing requires estimating the 
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matrix elements of the various operators in C e ^ between initial and final state 
mesons. The estimation of such hadronic matrix elements is notoriously difficult, 
and is generally accompanied by large uncertainties due to long-distance, non- 
perturbative strong-interaction physics. There are two factors which mitigate the 
effect of these hadronic uncertainties in the current phenomenological context: 

1) The dominant supersymmetric contribution to K° — K° and B° — B° mixing 
in the Minimal Supersymmetric Standard Model (MSSM) 21 results from a 5m 2 LL 
mass insertion, which gives rise to the same (V — A) 2 operator structure found 
in the Standard Model. This makes comparison of supersymmetric and Standard 
Model contributions relatively straightforward; 

2) For the B° — B° system, the valence quark approximation employed below 
is expected to be a good approximation. This belief is supported by lattice monte- 
carlo estimates of the B° — B° matrix element "fudge factor" Bb which give 
B B ^ l. 22 

The conventional result for the standard (V — A) 2 AS = 2 operator is 

(K^da^Wd^PLSjWK*) = 2 -fW K B a K (5a) 

where fx — 165 MeV is the K-meson decay constant and B^ = 1 corresponds to 
the "vacuum insertion" result. Various estimates of this matrix element place B^ 
in the range of 0.3 — 1.0, 23 with a value B^ ~ 0.7 favored by lattice gauge results. 22 
Matrix elements for the other hadronic operators appearing in £^L 2 follow from 
current algebra: 



(^°|[^P iS a]M /3 PL^]|^°) = ^f-^-) fWx (56) 

1Z \Tfl s -+- lud J 

(^ |[^^][^PL S J|^ ) = -^f-^^) f 2 K m 2 K (5c) 

12 \m s + raj J 



(K^dvPLS^dpPRS^K ) 
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{K°\[a a P L s^P R s a ]\K l 



m K 



m s + m di 



fK m K 



(5e) 



Similar fudge factors B b K e should be associated with each of the matrix elements 
(5b-e) above. Note that the matrix elements (5b-e) are enhanced relative to (5a) 
by a factor f Tr VL ) ~ 10. This enhancement, which is significant for (5b and 
d), combined with large numerical coefficients for these non-standard chiral con- 
tributions to £^5=2 (1)> raise the prospect of interesting enhancements to AMx 
outside of the MSSM, for which we will find that the 5fh 2 LL mass insertion gives the 
dominant contribution. The corresponding results for the AB = 2 matrix elements 



follow from (5) by setting E>b — ( 
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Phenomenological Upper Bounds on 5m 2 -, and <5m| , 

to 11 d A s B d A b B 




f^ rn d A s B \ 


f Sm d A ,s B , \ 




{ Ml ) 


{ Ml ) 


\ Ml )\ Mf ) 


{a(3){a'(3') 


AM K 




AM B 


(LL) 2 or {RRf 


(0.10) 2 


(0.0080)2 


(0.27)2 


(LL)(RR) 


(0.0060) 2 


(0.00049) 2 


(0.073) 2 


{LRf or {RLf 


(0.0082) 2 


(0.00066) 2 


(0.082)2 


{LR){RL) 


(0.044) 2 


(0.0035)2 


(0.14)2 



TABLE I. All numbers in the table must be multiplied by the factor (M^/TeV) 2 . Numerical 
values assume y/x ~ M g /M q — 1. Stricter (weaker) results generally apply to x < 1 
{x > 1). Bounds derived from ex assume maximal CP violation. Bounds from AMb must 
be scaled by (160 MeV// B ). 



Phenomenological constraints on the various FCNC mass mixings 5m 2 



d asb 
eff 



appearing in C, e £g =2 can ^ e derived by inserting the matrix elements (5) into £^,5=2 
and requiring AM™ y < AM^ P for each of the chiral contributions to the K^ — Ks 
mass difference: 



AM K ~ 2Mi2 = (K c 



m K 



C eff I K { 



3.52 x 10 _15 GeV 



(6) 



Barring accidental cancellations between the various chiral contributions, this con- 
straint is conservative in that the short-distance, Standard Model charm-quark 
contribution is itself of order AM^ P , and of the correct sign. Setting = 0.7 
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(only) and taking x — 1 yields the constraints on 5m| due to AMk shown in 
Table I. 

The results of Table I are significant and warrant comment. All the constraints 
from AMk are useful constraints (i.e. 5m\ B /M~ < 1), even for heavy squark 
masses Mq 3> 1 TeV. Indeed the constraint on 5m 2 LL , which is actually the weakest 
of the constraints from AMk, provides useful information on the down-squark mass 
matrix for Mq~ g as large as 10 TeV. Above 10 TeV, the upper bound on 5m| /M| 
rises above unity, at which point the perturbative expansion in powers of 5m 2 /M~ 
breaks down, and no further useful information remains on the form of the squark 
mass matrix. 

The constraint for 5rh 2 RR is identical to that for 5m 2 LL . However, Table I shows 
that in the presence of both types of mass insertions, even stronger constraints 
apply. This is because the numerical coefficients (e.g., 504) in £^5=2 (-0 anc ^ the 
operator matrix element (5) oc [m^/(m s + mrf)] are both enhanced relative to the 
5m 2 LL RR contributions. This results in useful constraints (i.e. 8m 2 LL RR /M~ < 1) 
on the squark mass matrix for sparticle masses M^g greater than 100 TeV. For 
the same reasons, the bounds on Smi and 5m| are similarly stringent and 
furnish useful constraints for sparticle masses M^g up to 100 TeV. 

The most stringent bounds on follow from the CP impurity parameter 

tu- Table I shows the constraints on the Sml which result from setting e s ™ sy < 

A dASB to K 

e e ^ p and assuming maximal CP-violating phases for each of the chiral contributions. 
Note that, in many cases, CP violation in K° — K° is sensitive to supersymmetric 
contributions from sparticle masses M^g > 1000 TeV. 

The analysis of AMg is identical to AMk, with appropriate substitution of 
flavor indices. Although numerically less stringent than those from AMk, the 
constraints in Table I can be equally restrictive for many supersymmetric models 
(like the MSSM) which predict 5m^ &B > «Jm| . 20 

In the MSSM with universal soft supersymmetry breaking, FCNC mass inser- 
tions are generated by renormalization effects between the unification scale Mjj 
and mw- 2 5 The resulting low-energy mass insertions are predominantly left-left, 
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and exhibit the following approximate flavor dependence: 

=CLL(t0,U)[VW u V] ij (7) 

where \ u is the diagonalized charge 2/3 quark mass matrix and V is the KM 
matrix. The renormalization-group coefficient cll is flavor independent, and is 
plotted in Figure 2 as a function of £o = m o/ m i/2 an d Ca = ^/ m i/2> where mo, 
m^2, A are respectively the universal soft supersymmetry-breaking scalar mass, 
gaugino mass and trilinear coupling at fj, — Mjj. Using the experimental upper 
bound 24 V t * d Vts < 7.3 x 1CT 4 on the KM angles associated with the dominant 
A 2 contribution, we observe that the radiatively-generated s is too small to 
contribute significantly to AM^ or e#. 

Such small FCNC results are a relatively unique prediction of the MSSM. 
Most nontrivial extensions of the MSSM contain extra Yukawa couplings, and 
these generically lead to large FCNC's. 6 In supersymmetric Flipped SU(5), 25 for 
example, we find potentially large FCNC's generated above the GUT scale: 20 

NLJy = -^2 Wily K^w) [3m ° + a2] (8a) 

NUJy = v™U*j ( 8& ) 

[^LdJy = -Vijv = ^[ A eA|A d + \ d \&I]ij ln (]j^) ( 8c ) 

where A6 is an a-priori unknown Yukawa coupling associated with a see-saw neu- 
trino mass mechanism. 25 

From the phenomenological constraints in Table I, we can derive useful 
bounds on the unknown GUT Yukawa A6- The strongest bounds come from the 
(Smj^ L )(Srh 2 RR ) contributions to AM^, AMg, from which we obtain: 
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e 2 + 6 ( M~ q x 



(96) 



S^ + aVlTeVy 



where we have used the approximate low-energy relation M| ~ (£q + 6)m^ 
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Equation (9) provides important knowledge about the unknown matrix A6 
and the pattern of soft supersymmetry breaking. For instance, if supersymmetry 
breaking takes the form of either mo or A, then (9) requires certain elements of A6 
to be quite small — smaller than expected from superstring theories, which relate 
Yukawa couplings to gauge couplings A ~ g ~ 0.7 at /j, — Mp\. One may be forced 
to conclude that these couplings vanish at the tree level in such theories — or that 
soft supersymmetry breaking takes the form of m^- 

Note that because of (8b), the Flipped contribution to CP violation ex van- 
ishes. Flipped SU(5) can, however, contribute significantly to "direct" CP violation 
e', Kl — > 7r°ee and Kl — > 7r°z/z7, as shown in ref. 20. 

We conclude that FCNC processes, and K° — K° and B° — B° mixing in 
particular, provide a sensitive probe of physics beyond the MSSM. In the presence 
of sparticles Mq, Mg & 1000 TeV, these FCNC processes provide an experimental 
window on the nature of physics at extremely high energies which should clearly 
be exploited to the fullest extent possible. 

The work of one of us (J.S.H.) was supported in part by the National Science 
Foundation under grant no. PHY-91 18320. 
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